Light hadron production in B c — > J/ip + X decays 
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Decays of ground state _B c -meson B c — > J/ip + mr are considered. Using existing parametrizations 
for B c — > J/ip form- factors and W* — > nil spectral functions we calculate branching fractions and 
transferred momentum distributions of B c — > J/ip + nn decays for n = 1,2,3,4. Inclusive decays 
B c J/ip + ud and polarization asymmetries of final charmonium are also investigated. Presented 
in our article results can be used to study form-factors of B c — > J/ip transitions, 7r-meson system 
spectral functions and give the opportunity to check the factorization theorem. 
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I. INTRODUCTION 



= 0.282 ± 0.038 ±0.074 



= 0.249 ± 0.045 +0 - 107 



Recent measurements of S c -meson mass and lifetime in CDF [l| and DO [|J experiments allow us to hope that 
more detailed investigation of this particle on LHC collider, where about 10 10 S c -events per year are expected, would 
clarify mechanisms of B c production and decay modes. Currently only products of _B c -meson production cross section 
and branching fractions of decays B c — > J/ipn, J/ip£v are known experimentally. For example, the following ratios 
r- — i are measured [jj: 
■ 

^! ob Br (B c ^ J/ije + u e ) 

On' (T B Br (B e -)■ J/ipK) 

i-C ■ for positron in the final state and 

! o Ba Br {Be -» J/j'^ + v l _ l ) 

>' <j B Br(B c ^ J/t/jK) -0.076 

on : 

for muon. These ratios are about an order of magnitude higher than the theoretical predictions based on current 
estimates of B c -meson production cross section and branching fraction Br(B c — > J/iplv) w 2% Q. The mode 
B c — > J/ipn was used mainly to determine precisely _B c -meson mass. No information on production cross section, 
decay branching fraction, and even the product of these quantities was determined in this experiment. 

Investigation of other jB c -meson decay channels and determination of their branching fractions will be one of 
interesting tasks of future experiments on LHC. Weak B c decays can be caused by decays of both constituent quarks. 
Dominant are c-quark decay modes, which amount to ~ 70% of all B c -meson decays. Unfortunately, none of such 
^ . reactions were observed, although large branching fractions are expected for some of these decay modes (for example, 
for B c —> B s p we have approximately 16% branching fraction). Mentioned above decays B c —> J/ipiv and B c — > J/ipir 
are examples of other class, caused by 6-quark decay. Total branching fraction of this process is about 20%. 

In the present paper we will fill the gap in existing theoretical predictions of _B c -meson decay branching fractions 
0-01 and consider multi-particle processes B c — > J/ip + nn with n — 1,2,3,4. These reactions are caused by weak 
6-quark decay b — > cW* — > cud and clean analogy with similar r-lepton decays (r — > v T + rm) can be easily seen. This 
analogy allows us to use existing experimental data on r-lepton decays and give reliable predictions of B c —> J/ip + mr 
branching fractions. 

In the next section we give analytical expressions for distributions of B c — > J/ip + nix decays branching fractions 
over invariant mass of the light hadron system and study different asymmetries of final J/-0-meson polarization as 
a function of this kinematic variable. In section III we use existing experimental data on r-lepton decays calculate 
branching fractions of B c — > J/tp + nir decays for n = 1,2,3,4. In section IV inclusive reaction B c — > J/ipud is 
considered in connection with duality relation. Short results of our work are given in the final section. 
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II. ANALYTIC RESULTS 

-B c -meson decays into light hadrons with vector charmonium J/ip production are caused by 6-quark decay b — > 
W* — ¥ cud (see diagram shown in figO}. The effective lagrangian of the latter process reads 

n eS = ^=v cb v: d [c+(fi)o + + C-^)0-}, 

where Gf is Fermi coupling constant, Vij are the elements of CKM mixing matrix, C±(/x) are Wilson coefficients, 
that take into account higher QCD corrections and operators 0± are defined according to 

0± = (d i u J ) v _ A (c i b 1 ) v _ A ±(d j u i ) v _ A (c i b j ) v _ A . 

In this expression i,j are color indexes of quarks and {q~\q2)v-A = 9i7^(l — 75)92- Since in our decays light quark 
pair should be in color-singlet state, the amplitude of the considered here processes is proportional to 

01 (M) = ^[(JV c -1)C + (m) + (JV c -1)C_(m)] 

If QCD corrections are neglected, one should set aif/i) = 1. Leading logarithmic strong corrections lead to dependence 
of this coefficient on the renormalization scale \x [8|], and on \l ~ mj, it is equal to 

Oi (mb) — 1.17. 

The matrix element of the decay B c — > J/ip + 7Z, where 1Z is some set of light hadrons, has the form 

M [B c ->• W*J/ip ->• TIJ/ip] = ^E^La^ef. (1) 
In this expression e K is the effective polarization vector of virtual VF-boson and 



= {J/4>\&y ll (l-T B )b\B e )=V ll -A». 

Vector and axial currents are equal to 

V p = {J/4> \c lfi b\ B c ) = ie^ef (p + k) a qpF v (q 2 ) , 

A„ = ( J/V \H^b\ B c ) = ef t F Q A (q 2 ) + (e*p) (p + fc)„ (q 2 ) + (e*p) q,F A (q 2 




Figure 1: B c -> J/ip + TL 
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Table I: Parameters of B c -meson form-factors 



where p and k are the momenta of B c - and J/f/j-mesons, q = p — k is the momentum of virtual VF-boson, and Fv(q 2 ), 
F^ ± (q 2 ) are form-factors of B c — > decay. Due to vector current conservation and partial axial current 

conservation the contribution of the form- factor are suppressed by small factor ~ (m u + m^) 2 /AfJ, , so we will 
neglect it in the following. 

One can use different approaches when deriving the form of the form-factors F(q 2 ). First of all, it is clear, that 
quark velocity in heavy quarkonia is small in comparison with c, so one can describe heavy quarkonia in the terms 
of non-relativistic wave- functions. This fact was used on the so called Quark Models [1, I9l4l4j|. In the following 
we will refer to this set of form-factors as QM.The speed of the final charmonium in £? c -meson rest frame, on the 
other hand, is large, so one can expand the amplitude of the considered here process in the powers of small parameter 
Mji^p/Mb,., as it was done in papers [15l420| . In what follows, we will refer to this set of form-factors as LC. One can 
also use 3-point QCD sum rules to obtain the information on B c — > J/ipW* form- factors [1, [HI, HH, (SR). 

In our paper we use the following simple parametrization of form-factors 



where numerical values of parameters Fi(0) and M po i e are presented in table [I] 
The width of the B c — > J / iplZ decay is 

dT(B c ^J/^U) = l-^pialWn w e^d^(B c ^ J/^TZ) 
where Lorentz-invariant phase space is defined according to 

d 3 p t 



d$(Q-nn...p n ) = (27r) 4 5 4 (Q-^p 4 )n. 



2£;(27r) 3 ' 

It is well known, that the following recurrent expression holds for this phase space: 

da 2 

d<S>(B c ^ J/ipK) = -^-d<f>(B c ^>- J/i)W*)d$(W* -+U). 
Using this expression one can perform the integration over phase space of the final state 1Z: 



i J d$ (W* -> K) e^e^* = (q^ - g V) Pt (q 2 ) + Wb>Ph (?) . 

where spectral functions L (q 2 ) are universal and can be determined from theoretical and experimental analysis of 
some other processes, for example r — > v T lZ decay or electron-positron annihilation e + e~ — » 7Z. Due to vector current 
conservation and partial axial current conservation spectral function is negligible on almost whole kinematical 
region, so we will neglect is in our paper. Explicit expressions for spectral function for different final states 1Z are 
given in the next section. 
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Differential distributions of longitudinally and transversely polarized J/^-meson in B c — > J/tp + lZ decays can easily 
be obtained from presented above expressions. In the case of longitudinal polarization the polarization vector e*" is 
equal to 

*/i_n\_ M \o n n M 2 + M 2 - q 2 



v ' 2M V \ M 2 
where z-axes is chosen in the direction of J/tp movement, M and M v are B c - and J/^-meson masses and 



(M + M v ) - q 2 / (M - M v ) - q 2 



M 2 V M 2 



Differential distribution has the form 



dT [B c ->■ J/^a=o + K] _ G 2 F M 3 V 2 b a 2 (3 n , 2 , M 4 ( ( 2 4M 2 g 2 



+ 2/3 2 (Af 2 - M£ - g 2 ) i^F^ 4 } . 
In the case of transversely polarized vector meson has the form 



y/2' y/2' 

and the corresponding differential distribution is 

dT [B c -> J/^a=±i + _ G^Vj, 2 TC r , - 



'^Wpt (q 2 ) [\F A \ 2 + M^ 2 \F V \ 2 ± ^Re (F Q A F V )} 



dq 2 32 

It should be stressed, that the above expressions are universal and spectral function p!p (s) depends on the final state 
K. 

If the polarization if final vector meson is not observed, the ^-distribution is, obviously, 

dT [B c -» J/jj + K] = dT [B c -» J/V^a + ft] , s 

dg 2 dq 2 U 

y A=0,±1 y 

It is also useful to study some polarization asymmetries. For example, polarization degree a is defined according to 

_ dT x=+ i + dTx=-i ~ 2dr A ^p 

dT x=+1 + rfr A= -i + 2rfr A=0 ' 
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q 2 « q 2 « (M Bc - Afy) 2 




Figure 3: Kinematics of B c — > J/tpud decay 



Production of transversely polarized, longitudinally polarized and unpolarized J/^-meson corresponds to a = f, 
a = — 1 and a = respectively. We would like to note, that in the framework of factorization model this asymmetry 
does not depend on final state TZ. So, experimental investigation of this asymmetry can be used for determination of 
B c -meson form factors and test of QCD factorization. In figJ5]we show g 2 -dependence of this asymmetry for different 
sets of B c -meson form-factors. One can easily explain qualitatively the behavior of these curves. Let us consider q 2 - 
dependence of asymmetry a in B c — > J / tpud decay. At low q 2 the direction of u- and (i-quarks momenta in i? c -meson 
rest frame will be close to each other and opposite to the direction of the momentum of J/^-meson. The spin of light 
ii-antiquark (d-quark) is directed along (opposite to) its momentum (see figl^), so quark-antiquark pair has A = 
projection on Oz axis. From angular momentum conservation it follows, that J/?/>-meson should also be longitudinally 
polarized. This can be observed in figure [21 where at low q 2 we have a = — 1 for all sets of £? c -meson form- factors. In 
high <7 2 -region, on the contrary, direction of quark and antiquark momenta are opposite to each other and J/^-meson 
stay at rest in _B c -meson rest frame (see fig'Ob)- As a result, final J/V>-meson is unpolarized in this region and a — 0. 
Another example is transverse asymmetry 

dT x=1 - dT x =-i 
aT = dT ' 

This asymmetry also depends only on i? c -meson form-factors and its dependence on squared transferred momentum 
is shown in figHJ 

III. EXCLUSIVE DECAYS 

In this section we present differential widths and branching fractions of the decays B c — > J/ip + nir using presented 
above universal formula ([2]) and specific expressions for spectral function ( < ? 2 )- 




2 4 6 8 10 



Figure 4: Transverse polarization asymmetry qt of final J/^- meson in B c — > JW+K decays as a function of squared transferred 
momentum q 2 (in GeV 2 ). Solid, dashed and dot-dashed lines stand for SR 0,[H, QM [|, and LC [U respectively 
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A. B c — ¥ J/ipiv 



Let us first of all consider two-particle decays B c — > J/-07T and B c — > J/4P- 

In the case of B c —> J/i^tt decay the W* — > tt transition is expressed through leptonic constant f v : 



(tt \uj^ 5 d\ 0) = V^/^. (3) 

The numerical value of this constant can be determined from tt —> pv p decay width: f„ « 140 MeV. The spectral 
function, that corresponds to vertex ^ is 

ft (q 2 ) = 2f 2 5 (q 2 ) . 

Using this spectral function it is easy to obtain the following values of B c — > J/ipir decay branching fractions for 
different sets of form-factors: 

Br LC (B c J/iJjit) = 0.13%, 
Br QM (B c -> J Urn) = 0.17%, 
Br SR (B c -> J/ipir) = 0.17%. 

B. B c -> J/4 + 2tt 

The 2n channel is saturated mainly by B c — > J/ipP decay. The W* — ¥ p transition vertex is also expressed through 
p-meson leptonic constant 

(p\vrf„d\Q) = V2f p M p e„ 

where f p ss 150 MeV. If one neglects the width of p-meson, the corresponding spectral function has the form 



p T (g 2 ) = 2f 2 S (q 2 - m 2 ) . (4) 

The branching fractions of B c — > J/ipp for different sets of form-factors are: 

Btlc (B c -> J/4p) = 0.38%, 
Btqm (B c -> J/4p) = 0.44%, 
Br SR [B c -> J/4p) = 0.48%. 

In order to take p-meson width into account, one can use experimental data on r — > v T + 2ir decay. The differential 
branching ratio of this reaction is equal to 

dY(T^v T 1l) _ G 2 F (m 2 T ~q 2 ) 2 2 n 



dq 2 16mn T 



(m 2 T + 2q 2 )p^(q 2 ) 



This method was used by ALEPH collaboration to measure the spectral function p^iq 2 ) in the kinematically allowed 
region q 2 < m 2 [23| and can be approximated by the expression (see fig|5^) 



p 2 j? (s) re 1.35 x 10 



-3 



s-4m 2 \ 2 1 + 0.64s 



(s-0.57) 2 + 0.013' 



where s is measured in GeV 2 . In figiljD we show corresponding distributions dT (B c — > J/ 4 + 2ir) /dq 2 . Solid, dashed 
and dash-dotted lines in this figure correspond to form-factors SR, QM, and LC respectively. The branching fractions 
of the decay B c — > J/4 + 2tt are almost equal to B c — > J/4p decay branching fractions: 

Btlc (B c -> J/4^) = 0.35%, 
Bxqm {B c -> J/i/nnr) = 0.44%, 
Bt S R (B c -> J/ipmr) = 0.48%. 
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Figure 5: fig (a) — spectral function p^; (b) — V 1 dF(B c — > J/ip + 2ir)/dq 2 distribution for different sets of B c -meson 
form-factors. Solid, dashed and dot-dashed lines stand for SR 

0,[U, Q M 3. and LC HI respectively 




Figure 6: Spectral function and differential width for B c — > J/ip + 3-7T decay. Notations are the same as in figO 



C. B c J/ip + 3tt 



In the case of B c — > J/ip + 3tt decay (where 3ir stands for the sum of tt ir°ir° and 7r 



decay modes) the 



G-parity of the final state is negative. So we can expect, that this mode is saturated by axial- vector resonance a\. 
The width of this state is too large to neglect it, so we cannot use the expression similar to Q for W* — > 3tt transition. 
The corresponding spectral function can be determined from experimental and theoretical data on r — > v T + 3ir decay. 
In our article we use the following expression to approximate this function (see. fig|6ji): 



pft? (s) w 5.86 x 10" 



1 + 190.. 



(s- 1.06) 2 +0.48 



2 ' 



Distributions over q 2 for different sets of i? c -meson form factors are shown in fig|Bp. The branching fractions of 
B c —> J /if) + 3ir decay are 



Br LC (B e 
BrQM (B c 
Br SR (B c 



J/ip + 3tt) = 0.52%, 
J/ip + 3tt) = 0.64%, 
J/ip + 3tt) = 0.77%. 



D. B c -> J/ip + 4tt 

In the decay B c — > J/ip + 47r both 7r~7r 7r°7r° and ir~ir + ir~iT~ modes are possible in the following we consider 
the sum of these states. The kinematically allowed region in r — >• v T + Air decay is too small to determine the 
form of spectral function pjT, so it is more convenient to use energy dependence of Air production cross section in 
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Figure 7: Spectral function and differential width for B c ~ > J/ip + 4n decay. Notations are the same as in fig[5] 

electron-positron annihilation. It is easy to obtain the following expression for this cross section: 

Aira 2 



a (e + e -> 47r) 



pf{s). 



Spectral function calculated from experimental data [24| is shown in fig 17k and later we use the following 
parametrization: 



pf (s) w 1.8 x 10" 



s - 16m 2 \ 1 + 5.07s + 8.63s 2 



(s - 1.83) 2 + 0.61 



The distributions corresponding to this spectral function are shown in figUJ). The branching fraction for different sets 
of i? c -meson form-factors are 

Bv LC (B c ^ J/ip + An) = 0.26%, 
Br QAf (B c -> J/ip + An) = 0.33%, 
Bt S r (B c -> J/ip + 4tt) = 0.40%. 



IV. INCLUSIVE DECAYS AND DUALITY RELATION 



Let us know consider the inclusive decay B c — > B c + X where X stands for an arbitrary state of light hadrons. On 
quark level this reaction corresponds to B c — > J/ip + ud decay. If one neglects u- and d-quark masses, the spectral 
function of W* — ¥ ud transition is energy independent and equals to 



Pt 



1 



In fig[8] distributions of B c — > J/ip + ud decay branching fractions for different sets of B c - meson form- factors are 
shown. Integrated branching fractions of this decay are 

Bt lc {B c -> J/xf, + ud) = 7%, 
Bt qm (B e -> J/ip + ud) = 8.6%, 
Br SjR [B c -> J/ip + ud) = 12%. 

It should be noted, that sum of presented above branching fractions (that is Br(B c —> J/ip + mr), n = 1,...4) 
gives only about 30% the inclusive decay branching fraction. So one could expect noticeable events with multi-pion 
production in _B c -meson decays. 

Below ifX-production threshold only 7r-mesons can be produced in ud-pair hadronization, so the duality relation 
should be satisfied 



(2m K +A) 2 



1 dT(B c ->■ J/tpud) 



'Br(B c ->■ J/ip + mr) 
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Figure 8: Differential B c — > J/ipud branching fractions for different sets of _B c -meson form- factors. Notations are the same as 
in figE] 
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Table II: B c — > J/tpTZ decays branching fractions (in %) for different sets of B c -meson form-factors 

where A is the duality window. If we restrict ourselves to n < 4 in the right-hand side of this relation, it is valid for 

A fa 0.6 GeV. 

It is interesting to note that this value is almost independent on the choice of i3 c -meson form-factors and close to the 
value of duality parameter in gg — > J/ipcc and \b —> J/ipcc reactions [25l |26| . 

V. CONCLUSION 

In our paper we study exclusive and inclusive decays of Z? c -meson into light hadrons and vector charmonium J /if), 
that is the processes B c — > J/tp + ud and B c — > J/ip + nir where n = 1, 2, 3, 4. According to QCD factorization theorem 
the amplitude of these processes splits into two independent parts. The first factor describes the decay B c —> J/ipW* 
and one can use existing parametrizations of _B c -meson form-factors to calculate this amplitude. The second factor 
describes the fragmentation of virtual W^-boson. The information about these processes was taken from experimental 
distributions of multi-pion production in r-lepton decays and electron-positron annihilation. 

Our results are gathered in table [III where branching fractions of multi-pion production in B c — > J ftp + nir for 
different B c -meson form-factors are presented. The last column of this table contains the branching fraction of the 
inclusive decay B c —> J/i/j + ud. It is clear that up to i^if-production threshold only it- mesons could be produced in 
B c — > J/ijj + X decay, so some duality relation should hold. In our article it is shown, that to satisfy this relation it 
is sufficient to integrate the inclusive spectrum up to squared transferred momentum q 2 = (2rriK + A) 2 . It turns out, 
that A is almost independent on the choice of B c -meson form-factors and equals to ~ 0.6 GeV. 

The other interesting point are the polarization asymmetries of final J / ip- meson. In the framework of factorization 
model these asymmetries do not depend on the final state TZ, so one can use them to investigate form-factors of 
-B c -meson and to test the factorization theorem. In our paper we present the polarization degree a — (dTx/dq 2 — 
2dT L / dq 2 ) / {dT t / dq 2 + 2dT L /dq 2 ) and transverse polarization asymmetry q t = (dY\ = i/dq 2 — dl a=-i / dq 2 ) / '{dT / 'dq 2 ) 
for different sets of form-factors. 
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